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A comparison is made between the degree of spin polarization of electrons excited by one- and 
two-photon absorption of circularly polarized light in bulk zincblende semiconductors. Time- and 
polarization-resolved experiments in (OOl)-oriented GaAs reveal an initial degree of spin polarization 
of 49% for both one- and two-photon spin injection at wavelengths of 775 and 1550 nm, in agreement 
with theory. The macroscopic symmetry and microscopic theory for two-photon spin injection are 
reviewed, and the latter is generalized to account for spin-splitting of the bands. The degree of spin 
polarization of one- and two-photon optical orientation need not be equal, as shown by calculations 
of spectra for GaAs, InP, GaSb, InSb, and ZnSe using a 14 x 14 k ■ p Hamiltonian including remote 
band effects. By including the higher conduction bands in the calculation, cubic anisotropy and 
the role of allowed- allowed transitions can be investigated. The allowed-allowed transitions do not 
conserve angular momentum and can cause a high degree of spin polarization close to the band 
edge; a value of 78% is calculated in GaSb, but by varying the material parameters it could be as 
high as 100%. The selection rules for spin injection from allowed-allowed transitions are presented, 
and interband spin-orbit coupling is found to play an important role. 

PACS numbers: 72.25.Fe, 42.65.-k, 78.47.-|-p, 72.25.Rb 



I. INTRODUCTION 



The optical injection of spin-polarized electrons in 
semiconductors, familiar since the 1980s)i is again at- 
tracting attention, due in part to the potential utiliza- 
tion of a spin-polarized electrical current in a technol- 
ogy called "spintronics" :^'- It is well known that linear 
absorption of circularly polarized light in a semiconduc- 
tor produces spin-polarized electrons in the conduction 
band.^ This occurs as a result of the entanglement of elec- 
tron spin and motion caused by the spin-orbit coupling in 
the material; in the absence of spin-orbit coupling, there 
would be no net spin polarization of the excited carriers. 
For many common semiconductors, the highest valence 
states are in the degenerate heavy and light hole bands 
at the r point. Consequently, the highest degree of spin 
polarization that can be achieved is 50%. Such a situ- 
ation occurs when the photon energy exceeds the band 
gap, but is not large enough to excite carriers out of the 
split-off band. This can be understood from selection 
rules that result from the symmetry of the states at the 
P point i 

One way to increase the spin polarization of the in- 
jected electrons is to use materials where the degeneracy 
between heavy and light hole bands is removed by strain 
and/or quantum confinement, so that one can excite car- 
riers only from one band. From the symmetry of the 
states, one then expects 100% spin polarization. And 
indeed both theory^ and experiments'*'^ have shown a 
significant enhancement of the degree of spin polariza- 
tion. The spin polarization could also be increased by 
using compounds with crystal structures where there is 



no valence band degeneracy iSt^iS 

Spin injection can also arise from two-photon absorp- 
tion. For certain applications this may have advantages 
over one-photon spin injection due to a much longer ab- 
sorption depth, which allows spin excitation throughout 
the volume of a bulk sample. Two-photon spin injection 
has been investigated in lead chacogenides (PbTe, PbSe, 
and PbS), which are cubic, and have direct fundamental 
band gaps at the L points.^ High degrees of spin polar- 
ization in these materials have been predicted^ but not 
observediiSiiiiiSii^ Our focus in this paper is on semicon- 
ductors that have a direct fundamental band gap at the 
P point, such as GaAs. Based on arguments involving 
the conservation of angular momentum, it was recently 
suggested that 100% spin polarization could be achieved 
in unstrained bulk GaAs from two-photon absorption.— 
Earlier theoretical calculations, however, predict a two- 
photon spin polarization of no more than 64% for this 
class of cubic semiconductors iSiiSiiSiii 

In this paper, we report results of time-resolved pump- 
probe experiments that show the degrees of spin polar- 
ization from one- or two-photon absorption are in fact 
comparable for GaAs. We also discuss in detail the var- 
ious effects that can complicate the direct experimental 
comparison of the spin polarization obtained by one- and 
two-photon excitation. We present microscopic calcula- 
tions of two-photon spin injection that go beyond the 
spherical approximation made by earlier calculations. We 
show how the simple argument based on conservation of 
angular momentum breaks down, and examine the tran- 
sitions that give rise to the partial spin polarization. The 
calculated one- and two-photon degrees of spin polariza- 
tion are not equal for all materials, and we find that. 
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in fact, two-photon spin injection can be fully polarized, 
but only from transitions that do not conserve angular 
momentum. 

Optical transitions near the T point can be summarized 
with sketches such as those in Fig. ^ The symmetry of 
the states at the T point of a crystal with zincblende sym- 
metry is as follows. The conduction band (Fgc) is s-like 
with two degenerate spin states, while the valence bands 
are p-like. The p-like orbitals are coupled to the electron 
spin to form four states (the heavy and light hole bands, 
Fg^) with total angular momentum 3/2 and two states at 
lower energy (the split-off band, Fy^,) with total angular 
momentum 1/2 The levels corresponding to the split- 
off band are not shown in Fig. ^ The selection rules for 
the transitions between these states are the same as for 
the states of a spherically symmetric system:^ Thus they 
can be understood using angular momentum arguments. 
By applying the familiar selection rule that one-photon 
absorption of circularly polarized light with positive hc- 
licity (iT^) must change the projection of total angular 
momentum by one sees that only the two transitions 
shown in Fig.^ are allowed. An examination of Clebsch- 
Gordan coefficients reveals that the transition from the 
ruj = —3/2 state of the valence band to the nij = —1/2 
state of the conduction band is three times as probable 
as the transition from the rrij = — 1/2 state of the valence 
band to the rrij = +1/2 of the conduction band. Thus, 
near the band edge, one expects a value of 50% for the 
degree of electron spin polarization 



a) 



b) 



c) 



P 



(1) 



where 7V| {N^ ) is the concentration of electrons with 
spin down (up). 

The idea of angular momentum conservation was ap- 
plied to two-photon absorption by Matsuyama et alt^ 
They argue that because the total angular momentum 
of the two right circularly polarized photons is -1-2, only 
the transition from mj — —3/2 to rrij — -fl/2 is allowed 
(see also Fig. ^). Therefore, they suggest that even in 
a bulk semiconductor with degenerate valence bands the 
resulting electron spin polarization should be 100%, and 
indeed with an opposite sign with respect to one-photon 
spin injection. 

On the other hand, the degree of spin polarization due 
to two-photon spin injection has been calculated several 
timesSiiMMi using the 8 band Kane modelJ^ Ivchenko 
calculated the degree of spin polarization in the limit 
of large spin-orbit splitting. ^'^^ Arifzhanov and Ivchenko 
improved the calculation by allowing the split-off band 
to act as an intermediate state; they gave the degree of 
spin polarization at the band edge as a function of Eg/ A, 
where Eg is the band gap energy and A is the spin-orbit 
sphtting.^^ For GaAs, one infers a 51% degree of spin 
polarization from their results. Note that, in contrast 
to one-photon spin injection, the degree of spin polar- 
ization of two-photon spin injection near the band edge 
depends not only on the symmetry of the states, but 







1/2 +1/2 +3/2 



FIG. 1: Optical transitions in a bulk zincblende semiconduc- 
tor for circularly polarized light allowed by the selection 
rules: (a) for one-photon absorption, (b) for two-photon ab- 
sorption as suggested by Matsuyama et al.^ (c) two-photon 
allowed-forbidden transitions with a conduction band as an 
example of an intermediate state, and (d) two-photon allowed- 
allowed transitions for vanishing interband spin-orbit coupling 
and light incident along a (001) direction. The quantum num- 
ber rrij for the projection of total angular momentum on the 
light propagation direction of all states involved is indicated 
in the figures. The thickness of arrows and adjacent number 
in (a) and (c) express the relative transition probabilities. 



also on various material parameters. When only two- 
band transitions are included, a very simple expression 
for the two-photon degree of spin polarization has been 
given in terms of the conduction and valence band ef- 
fective masses, from which one infers for GaAs a spin 
polarization of 48%. 

In most III-V semiconductors, the next higher conduc- 
tion bands are p-like {Ttc and Fgc)!^^ The role that these 
higher bands play in two-photon spin injection has not 
previously been investigated. It is known that k • p mix- 
ing with these bands is responsible for cubic anisotropy 
of two-photon absorptioni2Si2i The higher conduction 
bands can also act as intermediate states in the two- 
photon amplitude. Such transitions are qualitatively dif- 
ferent than transitions within the set of bands nearest to 
the fundamental band gapi^SiSiSliS^ 

In section ^ we review the symmetry of two-photon 
spin injection and present our calculation including the 
higher conduction bands. In contrast to previous calcu- 
lations of two-photon spin injection^SiiSiiSiii our calcula- 
tion is not perturbative in fc. In section llTll we present the 
experimental comparison of one- and two-photon spin in- 
jection. In section ^3 we discuss the transitions respon- 
sible for the degree of spin polarization in two-photon 
absorption (Figs, ^c and d). In an appendix we derive 
expressions for the degree of two-photon spin injection 
due to so-called 'allowed-allowed' transitions. 
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II. CALCULATION OF TWO-PHOTON SPIN 
INJECTION 

For an electric field of the form E (t) = exp{—iujt) + 
c.c. (we sometimes write E^^ = E^e^), the two-photon 
spin injection rate can be written phenomenologically as 

gi ^ (^jklm^j^^k^U^rn* ^^^^^ ^ijklm jg ^ ^.g^^]^ 

pseudotensor symmetric on exchange of indices j and k, 
and on exchange of indices I and m; superscript lower- 
case letters denote Cartesian components and repeated 
indices are to be summed over»?^ For the point groups 
Td and Oh, appropriate to most cubic semiconductors, a 
general fifth rank pseudotensor has ten independent com- 
ponents. Forcing the j ^ k and I m symmetries, and 
the condition for reality of S, C'™^'' = (^^jfeim)*^ leaves 
three independent real components. 

We define C2A = -i^'"""' and (23 = ImC'""', where 
the indices a, h, and c denote components along the stan- 
dard cubic axes [100], [010], and [001]. Then the three 
independent real components are ReC2A, ImC2yi, and C,2B- 
In the standard cubic basis, the non-zero components of 
C, are 



^abbbc 


^abccc 


^bacaa 
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^bccac 
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= K2A 


^abcbb 


^accbc 
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^aabac 


^bbcab 


. ^cacbc 
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^babbc 


^cbcac 
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as well as those generated by exchanging j ^ k and/or 
I m, for a total of 48 components. 

The point group symmetry allows spin injection for 
linearly polarized light, associated with Im^2A- However, 
from a microscopic expression for C in the independent 
particle picture [see Eq. Q), lO, or lO below], one can 
show that C must be purely imaginary due to the time 
reversal properties of the Bloch states. One might expect 
deviations from the independent particle picture within 
an exciton binding energy of the band edge^^*^ In what 
follows, we assume the independent particle picture is 
valid, which leaves the two-photon spin injection specified 
in terms of two real parameters C2A and (^2B- 

The component of the spin injection rate along one of 
the cubic axes can be written compactly (with no sum- 
mation convention) as^ 



S' = 2i (E, 



e; 



y {C2A 



IE, 



+ (2C2B-C2a)|K 



(2) 

If the material were isotropic, the spin injection rate 
could be described by only one real parameter; (2A — 
2(^2B and the second term in Eq. Q would be zero. 

The cubic anisotropy means that the two-photon spin 
injection from circularly polarized light depends on the 
angle of incidence of the light relative to the cubic axes. 
For circularly polarized light incident along n specified 



by polar angles 6 and relative to the cubic axes, 

2C2S — C2A 



S ■ h = T2C2A\E^r 1 



4C: 



2A 



-/I 



(3) 



where / {9, 0) = sin^ (261) + sin^ {9) sin^ (2(/.). The upper 
(lower) sign is for right (left) circular polarization. The 
analogous equation for two-photon absorption is given by 
Hutchings and Wherrett.^^ Equation Q has extrema for 
light incident along (001) and (111) directions. Due to 
the cubic anisotropy, the net injected spin is not always 
parallel to n, although it is when n is along (001) or 
(111). In particular, for light along a (001) direction, 
jSj = '2C2a\E^\'^, while for light incident along a (111) 
direction, |S| = (4/3) (C2A + C2s) j^^^j^. 



A. Microscopic calculation 

We calculate the photoinjection rate of net electron 
spin density, S, using second order perturbation theory 
with the light treated classically in the long wavelength 
limit. We ignore interactions amongst the electrons, and 
between electrons and phonons. We take the photon en- 
ergy to be below the band gap, and twice the photon 
energy to be above the band gap. We neglect any spin 
polarization of the holes, since their spin relaxation times 
are typically very short. In the Fermi's Golden Rule 
(FGR) limit, the photoinjection rate is time-independent. 

Expressions for the two-photon spin injection rate un- 
der these assumptions have been given beforei2iiSiiL2& 
However, all previous calculations used semiconductor 
models in which all bands are doubly degenerate. In 
such a case, one finds that 



s = § E (ck|s|c'k)f](^);^f]^i=;„_^<5[2c.-c.,„(k)],, 

(4) 

where |nk) is a Bloch state with energy 7iti;„(k), is a 
normalization volume, S is the spin operator, ujnmi^) = 
LUn (k) — ujjn (k) , the prime on the summation indicates a 

f 2) 

restriction to pairs (c, c') for which ujcc' = 0, and 1 
is the two-photon amplitude 

0(2) _ f e \^s^ {Eu- Vc,» (k)) (E^ ■ v„,^ (k)) 
Uc^ J ^ c^™-c^(k) ' 



with v„.m (k) = (nk| v |mk) where v is the veloc- 
ity operator. It is well known, however, that in real 
crystals of zincblende symmetry the spin degeneracy 
is removed,^^''^'^ albeit with a small energy splitting. 
Since we are using a model that accounts for this spin- 
splitting,'^^ we must generalize the earlier microscopic ex- 
pressions. Such a generalization was recently discussed 
for one-photon spin injection^ 

If the spin-split bands are well separated, FGR gives 



c.f ,k 



-k| S jck) 



n 



(2) 
c, Ujk 



S [2lj - (k)] . (6) 
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However, in GaAs the splitting is at most a few meV for 
conduction states within 500 meV of the band edgei^ 
Since this is comparable to the broadening that one 
would calculate from the scattering time of the states 
(and also to the laser bandwidth for experiments with 
pulses shorter than 100 fs), spin-split pairs of bands 
should be treated as quasi-degenerate in FGR. Thus in 
place of Eq. Q or we use, 



27r 



X -{S [2lu - uj,,y (k)] -^5[2uj~ tJc'u (k)]} 



(7) 



where the prime on the summation indicates a restric- 
tion to pairs (c, c') for which either c' = c, or c and 
c' are a quasi-degenerate pair. The coherence between 
quasi-degenerate bands is optically excited and grows 
with their populations, as is the case with simpler band 
models that neglect spin splittingi24i^iiS4iL2& Using the 
time reversal properties of the Bloch functions, the ex- 
pression for C*-''''™ that follows from ^ can be simplified 
to give 



c 



ijklm 



4 27r 



X Im 



J2 E-^P^-^^w] 

c,c',v,k n.n' 

(ck| S' |c'k) (1/^'='" 



ylmjk^ /2 



{Unv (k) - LO) {lOu'v (k) - Uj) 



(8) 



where 

yjklm 



K,„, (k),v„,,„(k)}^'={v;„ (k),v;„ (k)} 



I Un 



and {vi,V2}'^ = {v\v{ + v{v\)/2. 

The photoinjection rate for the density of electron- hole 
pairs is 

^=^E|"elk''^[2^-^-(k)]. (9) 



L3 
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From Eqs. (TJ and 0, the degree of spin polarization, 
P, can be calculated, since 



P 



2 S • n 

n N 



(10) 



The sign of P is chosen so that a positive P corresponds 
to an excess of electrons with spin down, i.e. spin oppo- 
site the photon angular momentum. 

To evaluate the degree of spin polarization, we 
use a k • p model that diagonalizes the one-electron 
Hamiltonian (including spin-orbit coupling) within a 
basis set of 14 T point states, and includes impor- 
tant remote band effectsi^ Fourteen band models (also 
called five-level models) have been used to calculate 
bandstructureSf24i2£i2Si21 as well as lineas^Si^ and non- 
lineai"2ii2M£ optical properties of GaAs and other semi- 
conductors. Winkler has given a recent review of 14 band 



TABLE I: Model parameters. 





GaAs 


InP 


GaSb 


InSb 


ZnSe 


Eg (eV) 


1.519 


1.424 


0.813 


0.235 


2.820 


Ao (eV) 


0.341 


0.108 


0.75 


0.803 


0.403 


E'o (eV) 


4.488 


4.6 


3.3 


3.39 


7.330 


Ao (eV) 


0.171 


0.50 


0.33 


0.39 


0.090 


A" (eV) 


-0.061 


0.22 


-0.28 


-0.244 


-0.238 


Pq (eVA) 


10.30 


8.65 


9.50 


9.51 


10.628 


Q (eVA) 


7.70 


7.24 


8.12 


8.22 


9.845 


P,5 (eVA) 


3.00 


4.30 


3.33 


3.17 


9.165 


711- 


7.797 


5.05 


13.2 


40.1 


4.30 


72 L 


2.458 


1.6 


4.4 


18.1 


1.14 


73i 


3.299 


1.73 


5.7 


19.2 


1.84 


F 


-1.055 














Ck (meVA) 


-3.4 


-14 


0.43 


-9.2 


-14 



models4i The 14 states (counting one for each spin) com- 
prise six valence band states (the split-off, heavy, and 
light hole bands), and eight conduction band states (the 
two lowest, which are s-like, and the six next-lowest, 
which are p-like). The states are given in more detail 
in Appendix and except for the split-off hole states, 
they are shown in Fig. Qi. 

The model contains 13 parameters chosen to fit low- 
temperature experimental data. Of the two parameter 
sets discussed by Pfeffer and Zawadzki for GaAs, we use 
the one corresponding to a = 0.085 that they find gives 
better resultsi^ For InP, GaSb, and InSb, we use parame- 
ters from Cardona, Christensen and Fasal,^ The param- 
eters are listed in Tableland the notation is described 
in Appendix 1X1 For cubic ZnSe, we use the parameters 
given by Mayer and Rossler~ and a calculated value of 
C/sfS; we use = — 0.238 eV to give a conduction 
band spin-splitting that matches the ab initio calculation 
of Cardona, Christensen and Fasal.^^ There is more un- 
certainty in the parameters for ZnSe than in those for 
the other materials, but we include it as an example of 
a semiconductor with a larger band gap. 

Note that although remote band terms are included 
in the 14 x 14 Hamiltonian, we have neglected the re- 
mote band contributions to the velocity operator. The 
effect of these contributions on one-photon absorption 
was discussed by Enders et alm^ Removing the remote 
band terms from our Hamiltonian changes P for GaAs 
by at most 2%. Thus we feel justified in our neglect 
of the remote band contributions to the velocity opera- 
tor. Another contribution to the velocity operator v, the 
anomalous velocity term, h{(T x 'VV) / (Am^c^^ should 
be included if /c-dependent spin-orbit coupling is included 
in the 14 x 14 Hamiltonian. For the results reported 
in the following section, we have neglected fc-dependent 
spin-orbit coupling. To test whether this neglect is justi- 
fied, we have repeated the calculation for GaAs including 
such coupling only between valence and lowest conduc- 
tion bands and the associated anomalous velocity; the 
coupling is parameterized by Co = 0.16eVA^ (note that 
Co is distinct from the /c-linear term Ck)- It decreases 
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Excess photon energy (meV) 



500 



FIG. 2: Calculated degree of electron spin polarization P in 
GaAs. The solid (dashed) line is for two-photon excitation 
with light incident along a (001) ((111)) direction and the 
dotted hne is for one-photon excitation. The dash-dotted line 
is for two- photon excitation calculated with an 8 band (Ft^, 
Fgu, and Fee) spherical model. The inset shows P close to the 
band edge. The sign of P is given by Eq. llTOl . 



the two-photon P by sa 2% for excess energies between 
0.1 and 200 meV. The decrease increases for larger excess 
energy, reaching « 5% for an excess energy of 500 meV. 

Our two-photon spin injection calculation is similar to 
the two-photon absorption calculation of Hutchings and 
Wherrett.^"'^ We can reproduce their results by removing 
remote band effects, which they did not include. 



B. Calculation results 

The calculated degrees of electron spin polarization, P, 
are shown for GaAs, InP, GaSb, InSb, and cubic ZnSe in 
Figs.|2Has a function of excess photon energy, 2huj — Eg, 
where Eg is the fundamental band gap. We also show, 
for comparison, the degree of electron spin polarization 
due to one-photon absorptioni^ For each semiconductor, 
the one-photon degree of spin polarization is 50% at the 
band edge as expected from the F point selection rules. 

In GaAs, so long as the excess photon energy is less 
than the split-off energy (341 meV) and greater than 
about 50meV, there is a near equality of one- and two- 
photon P's. 

Close to the band edge however, there is a feature of 
the two-photon P that has not previously been identified; 
it is seen more clearly in the insets of Figs. |2}E1 The 
values of the two-photon P at the band edge for each 
material are listed in Table ^ We discuss this feature 
further in Sec. IIVI but we note here that it does not 
appear in a spherical approximation. To show this, we 
have calculated the two-photon P with the 8x8 Kane 
model that includes only the valence bands and the Fgc 
conduction bands, and with 73 set equal to 72 and = 
0; the result, which is independent of crystal orientation. 



a, 




500 



500 



0.01 0.1 1 10 100 

100 200 300 400 500 

Excess photon energy (meV) 

FIG. 3: As Fig.El but for InP, GaSb, and InSb. 



TABLE II: Calculated band-edge two-photon P. 





GaAs 


InP 


GaSb 


InSb 


ZnSe 


[001] 

[111] 


-20.5% 
-60.0% 


58.7% 
-16.6% 


-58.9% 
-78.4% 


-49.3% 
-73.3% 


-14.5% 
-57.1% 



is shown in the dash-dotted line in Figs. EEl 

Both one- and two-photon P's decrease as the excess 
photon energy is increased. This is due to band mixing 
away from the F point, which changes the selection rules. 
At excess photon energies above the split-off energy, the 
one-photon P decreases due to transitions from the split- 
off valence band.ii The two-photon P also decreases due 
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FIG. 4: As Fig. H but for cubic ZnSe. 



to these transitions, but less so. 

The possibility of cubic anisotropy in two-photon spin 
injection was first pointed out by Ivchenko,^ although it 
has not been calculated until now. Cubic anisotropy in 
two-photon absorption, on the other hand, has been cal- 
culated by Hutchings and WherrettpSi They found that 
near the band edge two-photon absorption of circularly 
polarized light in GaAs should be about 10% greater for 
light incident along [111] compared to along [001]^ The 
results of our calculation for GaAs indicate that two- 
photon spin injection varies with crystal orientation by 
a similar amount. Hence the degree of spin polarization 
shown in Fig. [3 which is the ratio of the two, varies with 
crystal orientation by only a few percent for most pho- 
ton energies in the range we investigated. This is not the 
case, however, for excess photon energies very close to 
the band gap, as can be seen in the inset of Fig. |21 The 
cubic anisotropy is more substantial for ZnSe and InP. 



III. EXPERIMENTAL COMPARISON 

To experimentally measure the degree of electron spin 
polarization we performed a polarization-resolved pump- 
probe experiment, where the transmission of the probe 
pulses is measured as a function of the delay between 
w 150 fs circularly polarized pump and probe pulses. 
Specifically, we measure the differential transmission 
AT/T = {Je - Tq)/Tq, where Te (Tq) is the transmis- 
sion with (without) the pump. If the absorbance change 
induced by the pump is small (i.e., AaZ <C 1, where I is 
the sample thickness and Aa = aE — Oi^ is the differ- 
ence between the absorption coefficient with and with- 
out the pump, respectively), the differential transmission 
will be proportional to — AaL Furthermore, if this weak 
absorption change is caused by phase-space filling asso- 
ciated with a thermalized nondegenerate distribution of 
carriers, then the differential transmission will be propor- 
tional to the carrier density (AT/T cx N). These con- 
ditions are usually satisfied after about 0.5 ps for thin 



samples, low carrier densities and relatively high tem- 
peratures. Finally, if the holes do not contribute signifi- 
cantly to the phase-space filling, then the degree of polar- 
ization can be experimentally determined by measuring 
the differential transmission for pump and probe pulses 
having the same {AT/T)++ and opposite (Ar/T)+" cir- 
cular polarizations. For probe pulses near the band edge, 
{AT/T)++ cx 3Ni + N^ and {AT/T)+~ cx SN^+N^, as a 
result of the same selection rules described above in Sec. 
P Defining, 



^JAT/T) 



- [at/tY 



{AT/TY^ + {AT/Ty 



(11) 



we then have Pcxp — P- However, if these restricted 
conditions are not met, then Pexp may not directly yield 
the degree of polarization. The degree of polarization 
and the spin relaxation time may still be extracted, but 
other effects may have to be considered.*'*" 

Polarization-resolved differential transmission mea- 
surements were performed using pulses from an optical 
parametric amplifier (OPA)— pumped by a regenera- 
tively amplified Ti:sapphire laser operating at 250 kHz. 
The laser system was tuned to produce ^ 150 fs pulses 
at 1550 nm (signal) and 1650 nm (idler). Two beta bar- 
ium borate (BBO) crystals were used to generate 775 nm 
pulses from the signal beam and 825 nm pulses from 
the idler beam. The second-harmonic and fundamental 
pulses were then separated using dichroic beamsplitters. 
Thus, we used 775 nm pulses to excite the sample by one- 
photon absorption, 1550 nm pulses to excite the sample 
by two-photon absorption and 825 nm pulses to probe the 
transmission of the sample. 

We used a semi-insulating (impurity level less than 
10*^ cm~'^), 1 /j,m thick sample of [001]-grown bulk GaAs 
that was van der Waals bonded to the glass substrate. 
The experiments were performed at a temperature of 
80 K. Consequently, the probe beam was resonant with 
the band gap energy Eg, and the pump beams had 
an excess energy {2fiLo — Eg) of 90mcV, which is con- 
siderably less than the spin-orbit splitting energy of 
341 meV. The peak irradiances of pump pulses were 
~ 2.3GW/cm^ (fluence ~ 320mJ/cm-^) for two-photon 
and w llMW/cm^ (fluence w 1.2mJ/cm^) for one- 
photon absorption, exciting in both cases a carrier den- 
sity of « 6 X 10*^ cm""^. Probe pulses were 10 times 
weaker than pump pulses for one-photon excitation. 

The design and implementation of the experiments in- 
volved a substantial effort to remove all possible exper- 
imental artifacts that could influence direct comparison 
of the results obtained for optical pumping by one- and 
two-photon absorption. First, by using a pump-probe 
technique, we directly measured the degree of spin po- 
larization of electrons at fixed times after the generation 
process. Thus our data are more credible than that of 
experiments using time-integrated methods, where the 
measured degree of spin polarization has to be corrected 
using the ratio of the lifetime and the spin relaxation 
time of electrons J**'!- Second, we used a relatively thin 
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FIG. 5: The dynamics of differential transmission AT/T af- 
ter excitation by circularly polarized pump pulses with the 
excess energy of 90 meV as measured using probe pulses with 
the same (upper curve) and opposite (lower curve) circular 
polarization for one-photon (a) and two-photon (b) excita- 
tion. 
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FIG. 6: The dynamics of the degree of spin polarization 
of electrons Poxp after one-photon (solid squares) and two- 
photon (open circles) excitation as computed from data shown 
in Fig. El 



layer of GaAs. One- and two-photon absorption have 
substantially different excitation depth profiles, the for- 
mer being considerably steeper than the latter. In or- 
der to directly use Eq. ((TT|l to determine the degree of 
spin polarization, it is necessary to keep the sample thin 
enough to ensure that Aal <C 1 for each process; how- 
ever, the sample must be kept as thick as possible to 
maximize the magnitude of AT/T. We selected a 1 /im 
thickness as a compromise between these two tenden- 
cies. For the sample temperature and the wavelengths 
used here, the transmissions T of pump pulses for the 
two-photon excitation were larger than 99%. For pump 
pulses for one-photon excitation T « 30%, and for the 
probe pulses T « 70% M: Third, by using a probe wave- 
length different from the pump, we were able to improve 
signal to noise by using a spectral filter to eliminate scat- 
tered pump light. Fourth, we precisely characterized the 
quality of the circular polarization of all optical beams 
used. Due to the finite spectral bandwidth of the fem- 
tosecond optical pulses and the quality of the quarter 
wave (A/4) plates, the optical beams were not 100% cir- 
cularly polarized, but instead consisted of both and 
(7". However, the quality of the circular polarization of 
both pumps beams was nearly the same. A nominal cr+ 
polarization state of pump pulses for one-photon excita- 
tion was 95% tr"*" and 5% <7~, while for two-photon exci- 
tation a nominal polarization state was 94% <7+ and 
6% a~ . As another check, we used a broadband quarter 
wave plate to monitor the helicity of all beams. Finally, 
to avoid problems due to possible sample inhomogeneity, 
we focused all three beams (the pump beams for one- 
and two-photon excitation, and the probe beam) to the 
same position on the sample using a single achromatic 
lens. Their mutual spatial overlap was checked using a 
pinhole. 

The results of the pump-probe experiment for one- 



photon excitation by a pump are shown in Fig. 
The upper curve corresponds to probing with a probe, 
while the lower curve was measured with a cr^ probe. 
The difference between the different polarization condi- 
tions is caused by spin-dependent phase-space filling as 
described above. The resulting electron spin polariza- 
tion P as a function of time delay is shown in Fig. (full 
squares). The decay of P is due to the randomization of 
the initial spin polarization. From the data depicted in 
Fig. we infer a time constant of w 200 ps, which is con- 
ventionally considered as half of the spin relaxation time. 
The dominant spin relaxation mechanism is probably the 
precession about anisotropic internal magnetic fields (the 
D'yakonov-Perel' mechanism) 42*^2*^^ For two-photon ex- 
citation by a cr"*" pump, the AT/T signal was the same 
as in the case of one-photon excitation as illustrated in 
Fig. The resulting values of spin polarization P are 
shown in Fig. 1^1 (open circles) as a function of time de- 
lay. We can now directly compare the results obtained 
for one- and two-photon excitation. First, note that af- 
ter excitation by cr^ pump pulses P has the same sign 
in both cases, in a clear contrast to the predictions made 
by in Matsuyama et al}^ (see Fig. ^). Second, the ini- 
tial values of P for both one- and two-photon excitations 
are, within the experimental error (including the non- 
ideal polarization state of optical pulses used), the same 
and equal to the theoretical value of 49% expected for 
both at these photon energies. 



IV. DISCUSSION 

The prediction of a 100% degree of two-photon spin 
injection mentioned in Section uses arguments famil- 
iar from spherically symmetric systems. At first it might 
seem incorrect to even apply these in cubic systems. For 
the crystal Hamiltonian is not rotationally invariant and 



8 



thus does not conserve angular momentum: The lattice is 
viewed as fixed and able to provide any amount of torque. 
However, the deviation from spherical symmetry is small 
in many cases, and hence angular momentum arguments 
should have approximate validity. Stated more techni- 
cally, since Td is a subgroup of Oh, which is a subgroup 
of the full rotation group, the Hamiltonian can be written 
as the sum of spherical, cubic and tetrahedral parts with 
the latter two treated as perturbations4i*^2i^i*SiS The 
8 band Kane model (even including remote band effects 
but with = 73L and — 0) is spherically sym- 
metric and is often used to describe many properties. It 
has been used, in particular, for earlier calculations of 
one- and two-photon spin inject ion lii^iiSiiSiii In a spher- 
ical model however, the transitions depicted in Fig. 
do not occur. By examining the possible intermediate 
states [i.e. band n in Eq. (O], we can see which tran- 
sitions do occur, and understand the transfer of angular 
momentum. 



A. Allowed-forbidden transitions 

When the intermediate state is in the same band 
as either the initial or final state (a so-called "two- 
band transition" ) , one of the photons causes an in- 
traband transition. These two-band transitions dom- 
inate two-photon absorption in GaAs^^ and indeed 
in most semiconductors i^^iSiiS^ They are 'allowed- 
forbidden' transitions because the intraband transition, 
proportional to the velocity of electrons in the band, is 
zero at the F point. Consequently, it is not possible to 
derive the two-photon degree of spin polarization using 
the states at the F point as can be done for one-photon 
excitation (Fig.^) or other two-photon transitions (Fig. 
UJi). Instead, one must go away from the F point and 
sum over all k directions. With this caveat in mind, we 
nonetheless give a schematic illustration of a two-band 
transition in Fig. One should bear in mind that, 
away from the F point, one cannot in general associate 
states in the heavy hole band with — ±3/2 and states 
in the light hole band with = ±1/2, since this is only 
true for k || z. It is essentially due to this complication 
that the sum over directions of k gives a two-photon P 
that depends on the details of the bands. 

The slower decrease of the two-photon P compared 
to the one-photon P at excess photon energies greater 
than the split-off energy can be understood from a con- 
sideration of two-band transitions. The one-photon P 
decreases in this regime due to the selection rules involv- 
ing transitions from the split-off bandi The same selec- 
tion rules apply to the interband part of the two-band 
transition, but the intraband part of transitions from the 
split-off band is much weaker than the transitions from 
the heavy and light hole bands, since the latter excite 
to states higher in the conduction band that have higher 
velocity. 

There are also allowed-forbidden transitions of the 



three-band variety {hh-lh-c, hh-so-c, Ih—hh-c, and Ih— 
so-c); in these cases, the inter- valence band matrix ele- 
ments can connect states of opposite spin. Their effect on 
the two-photon spin polarization approximately cancels 
out in GaAs, as one can see by comparing a calculation 
that neglects themii with one that includes themii& 

Within a spherical model, allowed-forbidden transi- 
tions must conserve angular momentum; two-photon ab- 
sorption with circularly polarized light must transfer two 
units of angular momentum to each electron-hole pair 
that is created. In order to understand how this leads to 
an incomplete spin polarization, one should form eigen- 
states of angular momentum, even away from the F point. 
Such states can be formed in a spherical model with en- 
velope functions over an expansion of Bloch states. Any 
treatment of electron angular momentum must then take 
into account both the cell-periodic part and the envelope 
function part of the electron wavefunction. It is the lat- 
ter that is neglected in the argument of Matsuyamaii^ 
We plan to return to a more detailed discussion of this 
issue in a future publication. 

Yet even without that analysis it is clear that, in a 
simple two-band spherical model consisting of a single 
spin degenerate valence band and a single spin degenerate 
conduction band, the two units of angular momentum 
are divided equally between the two parts of the electron 
wavefunction. This can be inferred from the fact that the 
envelope function for the relative motion of the electron 
and hole has one unit of orbital angular momentum (i.e. 
it is a p wave).^'' A two-band spherical model can be 
mocked-up from an 8 band spherical model by setting 
the heavy and light hole band masses equaliS^i^ Doing 
so with the formula for the two-photon P given by Bhat 
and Sipefii one sees that in that case the two-photon P 
is 50% at the band edge. More generally, the maximum 
two-photon P in a spherical model is 64%.^^'^^ 



B. Allowed-allowed transitions 

Allowed- allowed transitions are those for which both 
matrix elements in the two-photon amplitude Q are non- 
zero at the F point. Allowed-allowed transitions have 
a different frequency dependence than allowed-forbidden 
transitions. Near 2huj ^ Eg the former varies as {2fiLLj — 
EgY^"^ while the latter as {2huj — EgY^"^ . Hence, allowed- 
allowed transitions can dominate allowed-forbidden tran- 
sitions in a frequency range close to the band edge. For 
GaAs, however this range is only lOmeV.^^'^^ As seen in 
the 14 band calculation shown in the inset of Fig. 12 the 
two-photon degree of spin polarization in this range can 
be very different from the rest of the spectrum. These 
transitions are necessarily due to lower symmetry parts 
of the Hamiltonian; in a system with true spherical sym- 
metry one could not have a two-photon transition from a 
p state to an s state, since two-photon transitions cannot 
connect states of opposite parity. 

The selection rules for allowed-allowed transitions are 



worked out in Appendix ^ Consider first the sim- 
ple approximation of vanishing interband spin-orbit cou- 
phng A~ (denoted A in Ref 31). Then the basis 
states given in Appendix ^ are the energy eigenstates 
at the r point. For (T+ polarized light incident along 
[001], the only allowed-allowed transitions are depicted 
in Fig. these can be derived from Table III of Lee 
and Fan.^^ The product of the two matrix elements in 
the two-photon amplitude is the same for both tran- 
sitions. Thus, if the spin-orbit splitting of the up- 
per conduction bands Ag can be neglected compared 
to the other energy differences, then P is zero [see Eq. 
(|B7ll with A~ = 0]. For cr+ polarized light incident 
along [111], the non-zero transitions are i) [Fg^, -1-1/2) to 
|r^„-3/2) to |r^„-l/2); ii) |f;,„,+3/2) to (|r^„-i/2) 
and |r'7„-l/2)) to |r^„+l/2); and iii) |r^„-3/2) to 
(|r;,„-l/2) and |F'7„-l/2)) to |F^„+l/2>. Here the 
prime indicates that the states are rotated so that the 
quantization axis is [111] rather than [001]. If the spin- 
orbit splitting of the upper conduction bands Aq can be 
neglected compared to the other energy differences, then 
the third of these is zero and the probability for the sec- 
ond is three times that of the first, resulting in P = —0.5 
[see Eq. (|B8^ ]. 

However, close to the band edge, where allowed- 
allowed transitions dominate, the full 14 band calcula- 
tion (see Table m or the insets of Figs. |21EJ does not 
agree with these simple arguments. There is significant 
difference between materials; for GaAs P = —0.21 and 
P = —0.60 for light incident along [001] and [111] re- 
spectively. The disagreement is due to the importance 
of the spin-orbit mixing between the valence and upper 
conduction bands, characterized by a nonvanishing A~. 

The interband spin-orbit coupling A~ would be zero if 
the material had inversion symmetry iSiS*^ In contrast 
to most of the other parameters in the 14 band model, 
the value of A^ has not been determined by directly 
fitting it to one or more experimental results. Rather, 
it has been calculated by various methods: the empiri- 
cal pseudopotential method (— 61meV for GaAs) j^i'^S'S 
the tight binding method (— 85meV for GaAs)?32, the 
ab initio linear-muffin-tin-orbitals method (— llOmeV for 
GaAs),^'^ and by an indirect fitting with a 30 x 30 k • p 
Hamiltonian (-70meV for GaAs).33-56 

In light of the variation in calculated values of the in- 
terband spin-orbit coupling A~ , we have investigated the 
dependence of the band edge two-photon degree of spin 
polarization on A^ . The result, shown in Fig.Ej is rather 
dramatic. First, it shows that for small A~, P due to 
[001] incident light is proportional to A^, whereas P due 
to [111] incident light is less sensitive to A~. Second, it 
indicates that a 100% degree of spin polarization could 
indeed be possible due to two-photon absorption. But 
this possibility is not due to the transfer of angular mo- 
mentum from the light to the electrons. Since it results 
from allowed-allowed transitions that are only non-zero 
due to the lack of inversion symmetry and could only oc- 
cur for certain crystal orientations, we suggest that some 




A" (eV) 

FIG. 7: Sensitivity of the GaAs band-edge two-photon P to 
A~. The solid (dashed) line is for two- photon excitation with 
circularly polarized light incident along a (001) ((111)) direc- 
tion as calculated with the 14 band model. The dotted line 
is Eq. fSTt . 



of the angular momentum comes from the crystal lattice 
itself. 

The selection rules for allowed-allowed transitions in- 
cluding interband spin-orbit coupling are given for [001] 
incident light in Eqs. (jBlHB6|l and an expression for the 
resulting spin polarization is given in Eq. (|B7|I . It is 
worth noting that P is independent of the valence-upper 
conduction momentum matrix parameter Q. 

This allows us to see how the small spin-orbit mixing 
between valence and upper conduction bands can have 
an important effect on the band edge spin-polarization. 
Allowed-allowed transitions between the unmixed states 
(Fig. ^) are proportional to the small matrix element 
Pq, which would be zero if there were inversion symme- 
try. And since the intermediate state is in an upper con- 
duction band, the energy denominator of the two-photon 
transition amplitude is large, which further reduces the 
amplitude for these transitions. The interband spin-orbit 
mixing, proportional to A~ / Eq, is small, but it intro- 
duces allowed-allowed transitions with a valence band as 
an intermediate state. Then, instead of being propor- 
tional to Pq, the transition is proportional to Pq, and the 
energy denominator is smaller. This allows the condition 
CPq ~ DPqA^ to be met with fairly modest interband 
spin-orbit mixing. 

V. CONCLUSION 

We have shown experimentally that the degrees of spin 
polarization produced by one- and two-photon spin injec- 
tion are approximately equal in GaAs at an excess photon 
energy of 90meV. This was also recently confirmed ex- 
perimentally by Stevens et aZ.^^ where for (lll)-oriented 
GaAs the measured degree of spin polarization of elec- 
trons generated by one- and two-photon excitation was 
found to be the same, in accord with our measurements 
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for (OOl)-oriented GaAs. The experimental results agree 
with our theoretical calculations, and they are not at 
odds with angular momentum conservation. As well, we 
have calculated the degree of spin polarization in other 
materials to show that the one- and two-photon degrees 
of spin polarization need not be equal. 

We have presented the first calculation of two-photon 
spin injection that goes beyond a spherical model. The 
cubic anisotropy of the two-photon P is small for most 
of the semiconductors we investigated at photon energies 
where allowed-forbidden transitions dominate, although 
it is somewhat larger in ZnSe and InP than in the oth- 
ers. AUowed-allowed transitions, which do not appear in 
a spherical model, and hence do not conserve angular mo- 
mentum, are found to strongly modify the two-photon P 
close to the band edge, and cause a large cubic anisotropy. 
We have identified the selection rules responsible for these 
transitions, and found that interband spin-orbit coupling 
plays an important role. 

Measuring the two-photon P due to allowed-allowed 
transitions would be challenging in most semiconductors, 
since they only dominate in a narrow energy range, and 
the absorption rate is small close to the band edge. How- 
ever, such a measurement could serve as means of de- 
termining the parameter A~, which contributes to the 
electron g factor^'^ and the spin splitting of bandsi^iiSi 

One should bear in mind that Coulomb effects may 
modify the two-photon P within an cxciton binding en- 
ergy of the band edge. The electron spin lifetime will also 
be shorter within an exciton binding energy of the band 
edge due to the Bir-Aronov-Pikus mechanism of spin 
relaxatiouiiS However, the energy range where allowed- 
allowed transitions dominate could in fact be larger than 
the exciton binding energy in some materials. 



IPs,, ±1/2) = i\S)\a±) 

|Py„±l/2) = ±i=|Z')|a±)-f i=|X'±zF')K) 

IPs., ±1/2) = T\|l\Z')\ai) + ^\X'±^Y')\a^) 

IPs., ±3/2) = ±^\X'±iY')\a±), 

where under the point group T^, \S) transforms like Pi, 
while {\X) ,\Y) ,\Z)} and {|X') , |y') , 1^')} transform 
like P4.^ 

At the P point, the energy between Pgc and Pg^, bands 
is Eg, the energy between Pyc and Pg„ bands is E'q, 
the energy between T^v and Tj^ bands is Aq and the 
energy between Pgc and Pyc bands is Aq. The mo- 
mentum matrix elements are Pq = i(h/m) {X\p^ \S), 
P^ = i{h/m){X'\p'=\S), and Q = i{h/m) (X'lpV \Z), 
where m is the electron mass. The interband spin-orbit 
coupling is 



3ih 



(Z'l (VI/ X p)^ |X) , 



and its sign has been discussed by Cardona et al. The 
parameters 71L, 72L, and 73^ are the usual Luttinger 
parameters that account for remote band effects on the 
valence bands. Since the 14 band model accounts for 
the Pgc, Ptc and Pgc bands exactly, modified Luttinger 
parameters are used in the 14 x 14 Hamiltoniani^ The 
parameter F accounts for remote band effects on the con- 
duction band (Pec), essentially fixing its effective mass to 
the experimentally observed value. Finally, the param- 
eter Cfc is the small fc-linear term in the valence bands 
due to interactions with remote bandsi^ 
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APPENDIX A: NOTATION 

The basis states for the 14 band model are (with 

|a+) = IT) and |a_) - U)), 

|P7„,±l/2) = ±^\Z)\a±) + ^^\X±iY)\a^) 

|P8„,±l/2) = T\l\\Z)\a±) + ^\X ±iY)\a^) 
1 



|Pg,,±3/2) = ±-^\X±iY)\a±) 



APPENDIX B: ALLOWED-ALLOWED 
CONTRIBUTION TO TWO-PHOTON SPIN 
INJECTION 

To calculate S and N due to allowed-allowed transi- 
tions, we can approximate all the matrix elements and 
energies in the two-photon amplitude by their value at 
the P point, thus avoiding the integral over k. Since the 
bands are doubly degenerate at the P point, we can use 
Eq. Q. 



1. Light incident along [001] 

Since we use a basis of states with spin quantized along 
z, (c, P| |c',P) oc (5c, c' and we have 



c'\,v,T 



(2) 

ci,j),r 



(k)] 



where c t and c J, are shorthand for the bands with states 

|r6c,±i/2). 
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For (7+ light, with polarization e^^ = (x + iy) /\/2, and 
S II z from Eq. (0) and the degree of spin polarization is 



P = 



y 






2 




2" 


y 




ci,ti,r 


2 

+ 




2" 



All but the Fgc states are not eigenstates at the F point 
due to spin-orbit coupling between upper conduction and 
valence bands parameterized by A~. The Hamiltonian 
at the F point in this basis has off-diagonal elements, 
but the order of the basis can be arranged so that it is 
block diagonal with blocks at most 2x2. For the bands 
|F7^,-hl/2) and |F7c,-M/2) (or for the bands |F7„,-l/2) 
and |F7c, -1/2)), the block is 



'Eg - Ao -2A-/3 
-2A-/3 E'q - Eg 



Since A"/ {Eq + Aq) <C 1, the off-diagonal part can be 
treated perturbatively. To first order in the perturbation, 
we have eigenvectors 



\soUl) = |F7.,±l/2) + 
\scUi) = |F7c,±l/2)- 



2A- 



1 



3 E'q + Ao 
2A- 1 
3 E'o + Ao 



|r7c,±l/2) 
|r7.,±l/2). 



For the Fa bands, the blocks are 



-E„ 



A-/3 



A-/3 E^-Eg + A'„ 
with eigenvectors to first order in A^/ {Eq + Aq) 

A- 1 



hhUi) - 


|F8,,±3/2) - 


\lhUi) = 


|r8.,±i/2)- 


hcUi) = 


|F8c,±3/2) + 


\lcUl) - 


|r8c,±l/2) + 



3 


E'o 


+ 


a; 


A- 




1 




3 


K 


+ 




A- 




1 




3 


K 


+ 




A- 




1 




3 


E'o 


+ 





-|F8c,±3/2) 



7ir8c,±i/2) 



-|F8„,±3/2) 



-|F8„±l/2) 



The non-zero matrix elements of • v in the eigenstate 
basis that can cause a two-photon transition between v 
and c are 



e^; • ^hci,ihi (F) 

e^j ■ ^djici (F) 
e,^ • VsciM^ (F) 

Guj ■ VcT.sci (F) 



^0 



A- 



3 E'q + A'q 



^0 + ^0 



(Bl) 

(B2) 
(B3) 
(B4) 



-A-Q/3 ( E'q + Aq , 
• ^soi..m (r) = ^..^ + 2 (B5) 



E'q+Aq \E'q + A'q 



/2 / 2A^ P' 



,(B6) 



where we have dropped terms second order in A~. Note 
that • ^hhi.ihi (r) = by an exact cancellation, as it 
should from symmetry considerations. Thus we have 



\\E^\^^^Q[APi + BPoA-] , 



where A = {E'q + A'q - Eg/2) ^ and B 
{E'„ + A'^y^ A/3. We also have 



n 



(2) 



where C={E'q^ Egl2)~^ 



\\E^\^ ^|^Q[CPl,- DPqA-] , 



D : 



1 



E'q + Aq 



Eg/2 + AQi \E'q + A[ 



and we have dropped terms proportional to QPq{A )^. 
The degree of spin polarization is then 

p ^ {AP^ + BPpA-f - {CP^ - DPpA-f 
{AP^ + BPoA-f + {CP/, - DPoA-f ' 



2. Light incident along [111] 

For (T+ light incident along [111], it is more tedious to 
obtain an expression like Eq. (|B7p since there are more 
non-zero matrix elements of e^^ • v than for cr"*" light in- 
cident along [001]. By rotating the basis to states quan- 
tized along [111], the matrix of elements of • v becomes 
simpler, but the Hamiltonian is no longer in 2 x 2 blocks. 
When A^ = 0, the latter is not an issue. In that case, 
we find 



(£) ^l^^ll^W 



(2) 



f-V 



' 3 



?i|£^^l';«QPo^(A + 2C) 
v3 



(2) 



'iQPiyilC-A) 



where A and C are as defined in the previous subsection, 
and t and J, are along [111]. With the assumption that 
A^ < - £'3/2, AfnC and we find that 



P(A- = 0) = -1/2. 



(B8) 
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